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Abstract. - We discuss the dynamics of a single hole in the t—J model at finite temperature, 
in the limit of large spatial dimensions. The problem is shown to yield a simple and physically 
transparent solution, that exemplifies the continuous thermal evolution of the underlying string 
picture from the T = string-pinned limit through to the paramagnetic phase. 



Following early pioneering work [1,2], and given added impetus by the discovery of high-X' c 
superconductivity, the study of single-particle excitations in magnetic insulators remains highly 
topical. At possibly its simplest, the dynamics of a single hole in an antiferromagnet (AF) is 
captured by the t — J model [3], itself the strong coupling limit of the Hubbard model and 
a paradigm for the physics of strong, local electron correlations: hole motion, with nearest 
neighbour (NN) hopping amplitude t, occurs in a restricted subspace of no doubly occupied 
sites, and in a background of AF coupled spins with NN exchange couplings J = At 2 /U . The 
string picture (see e.g. [4,5]), wherein the hole is pinned by the string of upturned spins its 
motion creates, has proven central in understanding the t—J model at T = 0: since spins in the 
resultant string are flipped relative to the Neel configuration, resulting for J > in energetically 
unfavourable exchange fields that generate a potential growing linearly with distance from its 
point of creation, the hole is thereby confined. And while string unwinding, and hence a finite 
hole mobility, is in general induced both by spin-flip interactions and Trugman loop motion 
[6], the underlying picture is robust even in d = 2 dimensions (see e.<7.[5]). 

The limit of large spatial dimensions, d — oo, has also proven highly fruitful in study of the 
t — J model [7-11]. Here the string picture has been shown to be exact, and the model solved 
at T = [7,8]; the additional roles of disorder [8], and of second NN hopping [9], have also 
been investigated. Moreover, in addition to providing exact solutions, the d°° limit serves as 
a starting point for l/d expansions, enabling a systematic approach to the fmite-d case. 

Given the success and utility of the large-d approach, we use it here to investigate another 
key element of hole dynamics in the t-J model: the role of temperature. That the spectrum 
of single-particle excitations, and hole conductivity, should change radically with T on the 
scale of the Neel temperature Tjy, is easily understood; for while the two-sublattice A/B 
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structure characteristic of the Neel ordered AF phase persists for < T < Tjv, the sublattice 
magnetization is not saturated: m(T) <ra(0) = l/2. In contrast to the T = limit, therefore, 
any site has a non-zero probability p a<J (where a= A or B distinguishes the sublattices) of being 
occupied by an electron of cither spin type, a = +/— or |/|. In consequence, single-step hole 
motion may incur either an exchange energy penalty — if a spin hops to the 'wrong' sublattice, 
as occurs inevitably at T = — or an energy gain if the converse occurs; and since spins in 
the string need not therefore be energetically penalized, a finite hole mobility results. Hole 
dynamics thus change character entirely on the scale of Tn itself, since in the paramagnetic 
phase for T>T N , where m(T) = 0, the energy barrier to hole motion vanishes, and the hole 
propagates essentially as a free particle. In this Letter we show that an exact, yet rather 
simple and physically transparent, solution to the problem can be given; providing in addition 
a benchmark against which approximations, such as the NCA [10,11], can be judged. 
The t— J Hamiltonian H = H t + H.j is given in standard notation by 

i,j 

the sums are over NN sites, and c\ a = cJ CT (l — hi- a ) embodies the constraint of no double 
occupancy by electrons. To ensure a well behaved limit for d = oo, the hopping and exchange 
couplings are scaled as 

t = U/y/Z , J=J*/Z (2) 

with Z — > oo the coordination number (Z = 2d for a hypercubic lattice). Spin- flip in- 
teractions are wholly suppressed for d°° (they contribute to 0(1/ d)), see eg. [7,8]; whence 
Hj—t^J^i j SfSj , simple molecular field theory (MFT) for which becomes trivially exact 
for d°° . To study hole dynamics we consider the thermally averaged, local one- hole Green 
function (retarded), given by 

G ii (w) = J2P(s)(i; S \(z-H')- 1 \i;s) ; (3) 

s 

\i; s) denotes the hole on site i with configuration s — {cfc} for the remaining spins, with 
energy E s under Hj\ H'=H—E S and z = uj + i0+. G,,(w) = (i;s\(z — if') _1 |i;s) is the local 
Green function for the given configuration s, the probability P(s) of generating which is given 
exactly for d°° by P(s) — Ilfc(^i) Pk(&k) (with corrections 0(l/d) ), reflecting the statistical 
independence inherent in d°° . Here p k (<r k ) , such that p k (a) = p aa for site k belonging to the 
a = A or B sublattice (with Pa^—Pb-^), is given by 

p aa = \ [1 + 2X a am(T)} (4) 

with X a = +l(— 1) for a = A(B); and the sublattice magnetization is given from simple MFT 

by 



m(T) = \ tanh 
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(5) 



with Tjv = j J* (fcs = l) the associated Neel temperature. 

To obtain Gu(u>) for an arbitrary configuration s, as now sketched, we first introduce a 
Feenberg self-energy Si(u) via 

G ii (co) = [z-S i (w)}- 1 . (6) 
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To obtain one may proceed thus: (a) Separate the resolvent operator G = (z — H') 1 as 

G = G Q + G HtG, where G a = (z-H' J )- 1 with H'j = H'-E s , such that (j;s'\G \i;s)=5 ij 5 ssl /z; 
(b) take matrix elements of G between (j;s'\ and |z; s) , obtaining thereby an 'equation of 
motion' for Gji(uj) = (j; s'\G\i; s) which (c) may then be iterated perturbatively in t. This 
generates the Nagaoka path formalism [1,2], enabling systematic construction of Si(u>) in 
unrenormalized form, i.e. as an explicit function of frequency lo. The Nagaoka path formalism 
is however a particular case of the Feenberg perturbation series (see e.g. [12]), the general power 
of which is to express Si(u) in renormalized form — i.e. as an explicit functional of the {Gjj} 
— enabling from Eq.(6) a self-consistent solution for the {Gjj}. 

The problem is particularly clear for the Z — > oo Bethe lattice, on which we focus explicitly 
(the hypercube can also be handled, the same physical ideas being involved but the algebra 
more complex). Here, Si(u>) is given by 

S l = Y / ^s\H t \j; S '}(j; S '\(z-HT 1 \j^')(j; S '\Ht\t; S ) , (7) 

j,s' 

where since hopping is solely NN, the spin configuration s', with energy E s > under Hj, differs 
from s only by a single hole/electron transfer; and {j; s'\H t \i; s) = —t. From the energy change 
involved in transferring a Cj-spin electron on site j in configuration s to an empty site i, it 
follows simply that E s > =E s +\j(TjU) p (T), reflecting either the energy cost (\j<jj=+l) or gain 
(XjO-j = —1) under single-step hole motion; here 

uj p (T) = J*m(T) (8) 

is the magnitude of the energy incurred in hole transfer which, involving necessarily sites on 
different sublattices, is equivalent to the exchange energy cost for an on-site spin-flip from 
<7 = +(— ) to — (+) on an A(B) sublattice site. From this it follows straightforwardly that for 
Z — > oo, (j; s'\(z — i?') _1 |j; s') =Gjj(u) — \j(TjW p (T)), whence as sought Eq.(7) yields Si as a 
functional of the {Gjj}: 

^^G^-Vj^) (9a) 

3 

(with sites j NN to i). The thermal average of Si, denoted by S (and naturally independent 
of the sublattice index i), follows directly noting likewise that Ga = G is independent of i: 

S = i^fwGV - \ a (JUJ p (T)) (9b) 



Finally, we recognize that G = [z — Si] 1 reduces simply to G = [z — S] 1 for Z — > oo (the 
corrections being 0(Z~ 1 )). Hence the desired result 

G(uj) =[z-tl (p Bl G(u - w P (T)) + Pm G(u + w P (T)))] (10) 

where ps CT (= Pa-<t) an d w p (T) are given from Eqs.(4,8) solely in terms of the sublattice 
magnetization m(T), whose T and dependence follows explicitly from the simple MFT 
equation, Eq.(5). 

Eq.(10) is exact and, with Eqs. (4,5,8), closed. Further, although we have merely sketched 
the derivation here its physical content is quite clear, involving simply the thermal probability 
(p Ba factors) with which a NN site to the initially created hole is occupied by a cr-spin electron, 
together with the associated energy cost/gain (the u> =F u p {T) shifts) incurred in NN hole 
transfer. For T = 0, where m(T) = | and hence p Ba = 8 al and oj p (0) = \ J*, Eq.(10) reduces 
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to G(co) = [z — tlG(ui — ij*)] -1 , the well known exact d°° result from the T — string 
picture [8] (and for finite-ci an approximation often used in practice [13]), yielding the familiar 
discrete spectrum characteristic of the confined/localized hole. In the paramagnetic phase by 
contrast, T >T N — ij*, where m(T) = = uo p (T) and p aa = \ V a, a, Eq.(lO) reduces to 
G(w) = [z — t;j;G(u;)] _1 , producing the semielliptic spectrum of full width 4i* characteristic of 
free hole motion in the d°° Bethe lattice. 

As a typical example, Fig.l shows for J*/i* = 0.4 the resultant hole spectrum D(to) = 
— 7r _1 Im G(lo) at four different T's; the spectral sum rule dwD(u)) = 1 is always sat- 
isfied. From Eq.(10), the thermal evolution of D(u) is controlled by that of the sublattice 
magnetization m(T), Eq.(5). At the lowest T shown, T = 0.32V, m(T) remains close to 
saturation: m(T)/m(0) ~ 0.997. Hence D(u) resembles closely its discrete T — counterpart: 
the peak positions coincide, with only minor thermal broadening. Increasing T to 0.5TV 
produces a modest reduction in m(T) to ~ 0.96m(0), but an appreciable effect upon D(w), 
with increased broadening and the incipient formation of a background continuum. The latter, 
a precursor of free-hole dynamics, grows with increasing T; see Fig. 1(c) for T = 0.7T/v where 
m(T) ~ 0.83m(0). By T = 0.9T N where m(T) ~ 0.53m(0), the semiellipse characteristic of 
free hole motion is clearly emerging, Fig. 1(d), showing also that the 'shortest strings' — those 
to the low-w side of D(uj) — are the last to be thermally eroded. Finally, in the paramagnet 
for T > = \ the resultant semielliptic D(uj) is independent of T for obvious physical 
reasons. 

Fig.l shows clearly that while the effective spectral width is on the scale of the hopping t*, 
the essential character of the spectrum evolves on the typically much smaller scale of Tjy = \ J* 
(= 0.1t» in Fig.l). Qualitatively similar behaviour occurs on varying J*/i*; the rule of thumb 
being that for any fixed T/Tjv (and hence constant m(T) and p aa , see Eqs.(4,5)), increasing 
J* amplifies the hole transfer shifts in Eq.(10) (since u) p (T) = J*m(T)), thus producing a 
somewhat 'colder' spectrum. 

The behaviour described above is reflected also in the hole conductivity <t(lu;T), itself of 
0(1/ d) in high dimensions but with cr(w;T)/t 2 (or cquivalently Tr<rxx(oj;T)) of order unity 
and a much studied quantity [7,8] . Due to the absence of vertex corrections for d°° , the (real) 
dynamical conductivity a(u>;T) is expressible in terms of the hole spectrum D(u>); and an 
exact expression for it may be deduced along the lines sketched above for D(ui). We quote 
only the result (which is independent of lattice type), namely 

(1-e-^ /."dwie-^DW IpBiD^i+w-Wpjr)) ip^D^+w+Wp^)] 
<7(lu;T) OC -oo — (11) 

where (5 = l/T and extraneous constants have been dropped for clarity. Temperature enters 
Eq.(ll) in two ways: via (a) explicit Boltzmann factors and (b) the T-dependence of D(u>), 
p Ba and u) p (T). In previous work [8] the latter have been neglected entirely and replaced by 
their T = counterparts with p Ba = S„i and lu(0) = \ J*, whereupon Eq.(ll) reduces to the 
result of [8]. This amounts to retaining only the T = Neel spin configuration in calculating 
a(uj;T), but is not a physically realistic approximation: since D(u>) for T = is discrete, the 
dynamical conductivity consists of a series of (5-functions whose weight alone depends on T 
[8], and the dc conductivity <Jo(T) = <j(ij = 0;T) vanishes for all T. While the introduction of 
disorder [8] broadens the (5-peaks and leads to non-trivial behaviour, the absence of a proper 
T-dependence limits somewhat the utility of this approximation. 

Eq.(ll) by contrast is exact at finite-T. To illustrate it, Fig. 2 gives the resultant dc 
conductivity <to(T)/cto(Tjv) vs. T/Tm for J^/i* = 0.4. As expected it is Tn — \J*, and 
not the hopping t* (= lOT^v in Fig. 2), that sets the scale for the T-dependence of the 
conductivity in the AF phase. er (T) is non-zero for any T > 0, reflecting the fact that 
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true hole confinement/localization is exclusively a T — phenomenon, as embodied also in 
the corresponding loss of the discrete T = hole spectrum at finite T (Fig.l); the progressive 
increase of cr (T) with T in the AF phase likewise reflects the enhanced thermal broadening 
of D(ui). At T = Tjv the conductivity acquires a sharp cusp, across which dao(T)/dT changes 
sign. In the paramagnetic phase, p aa — |, ui p — and (see Eq.(lO)) D(ui) is independent 
both of T and J*: Eq.(ll) for a{uj;T) thus reduces in practice to its well known J* = limit 
[7,8]. For T > TV the T-dependence of cr (T) is thus controlled by the sole remaining scale, 
f*, diminishing monotonically with increasing T and approaching oo(T) <~ 1/T behaviour for 
T/t* > 1 — 2 (which asymptote is evident from the j3 = 1/T prefactor in Eq.(ll) for uj = 0). 

In summary we have shown that in the large dimensional limit, the finite tempera- 
ture dynamics of a single hole in an AF, as described by the t — J model, admits to 
a simple closed solution that captures in a physically transparent fashion the continuous 
thermal evolution of the underlying string picture, from the T — string-pinned limit 
to the paramagnetic phase T > TV where there is no exchange energy barrier to hole 
propagation. Although for clarity we have focussed explicitly on the Bethe lattice, the 
analysis is readily extended both to lattices with non-retraceable paths and, relatedly, to 
the case of a second NN hopping, which even for a Bethe lattice involves non-retraceable 
loop paths and likewise admits to a closed solution for all T [14]. The additional effects 
of disorder, as hitherto considered for T = [8], are also easily encompassed; this, to- 
gether with an elaboration of the present work, will be described in a subsequent paper. 
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Fig. 2. - Static conductivity a (T)/a (T N ) vs. T/T N for J„/U = 0.4. Note that T N = U/W. 



